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METHOD  OF  CALCULATING  THE  FREEZING  RATE  OF  SINGLE-COLUMN 
WATER-PERMEABLE  SOILS 

[Metod  Rascheta  Skorosti  Zamorazhivaniya  Fi 1 1 truyushchi kh  Gruntov 
Odinochnoy  Kolonkoy] 

A.  I.  Pekhovich,  Candidate  in  Technical  Sciences 

The  B.  Ye.  Vedeneyev  All-Union  Scientific  Research  Institute 
of  Hydraulic  Engineering 


I.  Introduction 

As  is  well  known,  the  artificial  freezing  of  soils  is  accomplished 
by  means  of  implanting  freezing  columns  in  which  brine  is  circulated 
in  the  ground  at  certain  intervals.  The  brine  is  cooled  by  cooling 
devices.  Initially  individual  frozen-soil  cylinders  are  formed 
around  each  column,  and  these  columns  then  fuse,  thus  forming  a 
solid  body.  Depending  on  the  nature  of  the  construction  operations, 
the  columns  are  placed  in  very  different  ways:  in  a straight 

line  or  in  a circle,  in  a single  row  or  in  several  rows,  etc. 

Predicting  soil  freezing  consists  primarily  of  selecting  the 
type  and  power  of  the  cooling  devices  and  the  dimensions,  number 
and  Si fang^ment  of  the  freezing  columns.  In  design  it  is  also  necessary 
to  develop  a prodticcizr  eraph  of  the  soil  freezing  operations  which 
must  indicate  the  operating  intflie  af  ♦■he  freezing  devices.  It  is 
quite  obvious  that  to  solve  all  of  these  p?80ic<r.c  and  to  select 
the  most  economical  version  the  designer  must  be  able  to  deccr-'^e 
the  soil  freezing  rate. 

However,  the  existing  methods  of  calculating  the  freezing  rate, 
of  water-permeable  soils  are  insufficiently  precise  for  practice, 
which  is  primarily  due  to  the  excessively  rough  allowance  for  the  influx 
f heat  from  the  unfrozen  ground.  Usually  the  thermal  influx  value 
av',opted  in  the  calculations  is  far  below  the  actual  value.  This  is 
CO.  firmed  by  experimental  materials  on  the  freezing  of  water-permeable 
soi..,s  obtained  in  recent  years  by  thq  Ice-Heat  Laboratory  of  the 
B.  Yv  . Vedeneyev  All-Union  Scientific  Research  Institute  of  Hydraulic 
Engineering  and  data  from  field  observations  of  a frozen-soil  cofferdam 
in  the  construction  of  * i>>aroelectric  power  plant,  as  well  as  the 
works  o/  B.  V.  Proskur>akov  [1]. 

The  {.-■  cicion  with  which  freezing  rate  is  calculated  is  also  reduced 
because  of  jnortcomings  in  the  methods  used  to  determine  the  heat 
absorption  factor  of  the  columns.1 


^ee  the  author's  next  article  in  this  collection. 


r 


This  article  analyzes  the  method  for  calculating  the  freezing  rate 
of  homogeneous  water-permeable  soils  by  a single  column;  the  method 
has  been  developed  in  works  of  B.  V.  Proskuryakov  [1]  and  I.  A. 

Charnyy  [2] . 

Our  task  can  be  formulated  as  follows. 

A single  freezing  column  of  radius  rQ  is  implanted  perpendicularly 

in  a homogeneous,  infinite  water-permeable  layer  of  ground  (Figure  1). 

The  base  and  top  of  the  layer  are  two  parallel  planes  which  are  permeable 
to  water  and  heat.  The  length  of  the  column  L is  equal  to  the  thickness 
of  the  layer  h. 

At  a considerable  distance  from  the  column,  the  filtration  flow 
has  velocity  v and  temperature  d.  The  ground  freezes  at  temperature 
d^.  The  physical  constants  of  frozen  and  unfrozen  ground  are  given. 

The  column's  heat  absorption  is  a known  function  of  the  frozen  ground's 
thermal  resistance.  It  is  necessary  to  determine  how  the  dimensions 
of  the  frozen  ground  depend  on  time. 


?- 


Frozen  ground  cylinder 

Unfrozen  ground  / Freezing  column 


Figure  1.  Freezing  of  Water-Permeable  Soil  by  a Single  Column. 

In  solving  this  problem,  we  will  introduce  the  following  basic 
assumptions,  the  influence  of  which  will  be  examined  later: 

1)  the  thermal  flux  which  changes  the  temperature  of  the  frozen 
ground  may  be  taken  into  account  by  increasing  the  computed  value  of 
the  latent  ice  formation  heat; 


I 


2)  the  temperature  distribution  in  frozen  ground  at  any  moment  of 
time  corresponds  to  the  established  thermal  state; 

3)  the  frozen  ground  has  the  shape  of  a cylinder; 

4)  the  axis  of  the  frozen  ground  cylinder  coincides  with  the  axis 


of  the  freezing  column. 

f 


Of  the  other  less  significant  assumptions  which  we  will  introduce 
but  whose  influence  we  will  not  examine,  let  us  note  the  following: 

1)  the  thermal -physical  constants  of  the  frozen  and  unfrozen 
ground  are  not  a function  of  temperature; 

2)  all  of  the  ground  water  freezes  at  the  same  temperature; 

3)  the  coefficient  of  heat  transfer  from  the  surface  of  frozen 
ground  to  the  freezing  water  layer  is  infinitely  great;  i.e.,  we  will 
ignore  the  transfer  heat  resistance  from  the  water  to  the  frozen 
ground  cylinder. 

2.  Thermal  Balance  Equation  and  Its 
The  thermal  balance  equation  of 

form: 

of  - Qc  - Qfi 

where 

Q.  is  the  thermal  flux  from  the 

l 

given  off  at  the  freezing  boundary; 

Qc  is  the  thermal  flux  which  passes  from  the  frozen  ground  to  the 
freezing  column  (the  thermal  absorption  of  the  column); 

Q^.  is  the  thermal  flux  which  passes  from  the  filtration  flow  to  the 
frozen  ground  cylinder; 

Q is  the  thermal  flux  due  to  the  reduced  heat  content  in  the 
g 

frozen  ground. 

Because  of  the  above-mentioned  assumption  that  the  thermal  flux 
can  be  taken  into  account  by  simply  increasing  the  calculated  value 

of  the  latent  ice  formation  heat,  thermal  balance  equation  (1)  is 
written  in  the  following  form: 


Components 

frozen  ground  has  the  following 

- V (1) 

latent  ice  formation  heat  which  is 


Of 


(2) 


The  thermal  flux  from  the  latent  ice  formation  heat  is  equal  to 


Of 


(3) 


where 


c is  the  amount  of  heat  given  off  when  a unit  of  ground  volume 
freezes; 


L is  the  length  of  the  freezing  column; 
r is  the  radius  of  the  frozen  ground  cylinder; 
t is  time. 


B.  V.  Proskuryakov  [1]  has  demonstrated  that  the  thermal  flux  which 
passes  from  the  unfrozen  water-permeable  ground  can  be  represented  by 
an  equation  with  the  following  form: 

Qfl--  8 ]/  iiSft L ('■»-»„ );  (4) 


X^  is  the  thermal  conductivity  factor  of  unfrozen  ground; 
cw  is  the  thermal  capacity  of  water; 

Yw  is  the  specific  weight  of  water; 
v is  the  velocity  of  the  filtration  flow; 
d is  the  temperature  of  the  filtration  flow; 
i?q  is  the  freezing  temperature  of  the  ground. 

The  last  term  in  equation  (2),  the  column's  heat  absorption,  is 
in  each  particular  case  a completely  determined  function  of  the  radius 
of  the  frozen  ground  cylinder: 

Qc  = f(r).  (5] 

3.  Solving  the  Thermal  Balance  Equation 

By  substituting  expressions  (3),  (4)  and  (5)  into  heat  balance 
equation  (2)  we  will  have 


where 


~~L'r  ~dS  —fir)- A | r , 


A=H  |/  ({)  — ;>(,). 


By  separating  variables  r and  t and  by  integrating  within  the 
limits  of  t = to  i = t?  and  r = r^  to  r = r2>  we  find 


;i=2 

■ I ( r ) ,11  r 
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The  integral  of  the  right-hand  portion  of  equation  (8)  can  be  most 
simply  determined  by  graphic  or  tabular  integration. 

The  maximum  size  of  the  frozen  ground  cylinder  r = r is  determined 

' max 

from  the  condition  that  the  thermal  absorption  of  the  column  is  equal  to 
the  thermal  influx  from  the  filtration  flow: 

f(r)  — A \r7.  (9) 

Computation  equation  (8)  which  is  used  to  determine  the  velocity 
of  freezing  is  correct  with  any  form  of  the  dependency  = f(r)1. 

The  uncertain  form  of  the  latter  [dependency]  made  it  possible  to 
present  computation  equation  (8)  in  integral  form,  which  is  very 
inconvenient  for  practical  use.  This  fact  causes  us  to  search  for 
short-cuts. 

The  course  of  freezing  consists  of  two  periods.  During  the  first 
period  the  temperature  of  the  brine  in  the  column  drops;  during  the 
second  period  the  temperature  of  the  brine  is  constant. 

The  first  period  changes  to  the  second  when  a frozen  ground  cylinder 
cylinder  is  formed  which  has  a radius  determined  from  the  equation: 

fir)~  (10) 


where 

rQ  is  the  radius  of  the  freezing  column; 

\ 2 is  the  thermal  conductivity  factor  of  the  frozen  ground; 

9m  is  the  maximum  low  temperature  of  the  brine. 

Equation  (10)  is  obtained  by  jointly  solving  equation  (5)  with  the 
expression  for  the  heat  flow  in  a cylinder  which,  as  is  known  from  heat 
transfer  courses,  has  the  following  form: 


We  will  determine  the  freezing  time  separately  for  the  first  and 
second  periods  of  freezing. 

During  the  first  freezing  period  we  will  assume  that  the  heat 
absorption  of  the  column  is  constant  (averaged);  if  desired,  it  is 


. 


lThe  method  for  determining  this  dependency  is  given  in  another  article 
by  the  author  later  in  this  collection. 


possible  to  refine  the  calculation,  for  which  we  should  determine 
the  freezing  time  in  sections  while  assuming  that  each  section  has  its 
own  constant  column  heat  absorption  value. 

During  the  second  freezing  period  the  brine  temperature  is 
constant,  and  therefore  the  heat  absorption  of  the  column  will  be 
expressed  by  equation  (11)  during  this  time. 

It  must  be  noted  that  in  a number  of  cases  the  entire  process 
of  freezing  takes  place  during  the  first  period,  and  in  some  other 
cases  the  first  period  comprises  an  insignificant  portion  of  the 
total  freezing  time.  In  these  cases  it  is  obvious  that  the  freezing 
time  is  determined  by  calculating  one  period. 

Thus,  let  us  determine  the  calculated  dependency  for  establishing 
freezing  time  when  the  column's  heat  absorption  is  constant. 

Heat  balance  equation  (6)  acquires  the  following  form: 


from  which  we  have 


‘IvLzr  ' 


— = Qc  - A \ r, 


r 


,1 
2 r.L 


dr. 


(12) 


In  order  to  integrate  the  latter  differential  equation,  we  will 
introduce  a new  variable 

z --  Vr, 

and  stipulating  that 


and 


(13) 


(14) 


we  find 


from  which  we  have 


d-  = Bx -22s  -2.VZ---2.Y,  W~- 
V A i j 

By  integrating  this  differential  equation  within  the  limits  of 
t = to  x = r2  and  : = ij  to  z = z,  and  by  returning  to  the  initial 

variable  r,  we  finally  find  the  calculated  dependency  which  is  used  to 
determine  the  freezing  and  thawing  rates  of  the  frozen  ground  cylinder 
when  the  heat  absorption  of  the  column  is  constant: 


-6- 


(15) 


~n  = -B 


• I TV ~r‘ /+'v«  ,r2~r*'l  +2‘Mr-:  ~r>  I 


+ 2,V;in 


ri  — A’ 


* V. 


where  is  the  initial  radius  of  the  frozen  ground  cylinder  at  time  t^. 

From  formula  (15)  it  is  evident  that  the  freezing  rate  noted 
when  the  column's  heat  absorption  is  constant  is  a function  of  only 
two  parameters:  and  N^.  This  makes  it  possible  to  plot  graphs 

which  are  convenient  for  making  the  calculations  (Figures  2 and  3). 


Figure  2.  Graph  for  Determining  the  Rate  with  Which  a 
Frozen  Ground  Cylinder  is  Formed  when  the  Heat  Absorption 
of  the  Column  is  Constant. 


Assuming  that  the  left-hand  portion  of  equation  (12)  is  equal  to 
zero,  we  find  the  maximum  radius  of  the  frozen  ground  cylinder: 


r 

max 


(lb) 


Let  us  now  attempt  to  determine  the  calculated  dependency  of  the 
speed  with  which  the  frozen  ground  cylinder  is  formed  when  the  brine 
temperature  is  constant,  which  is  what  occurs  during  the  second  freezing 
period.  In  this  case  the  heat  absorption  of  the  column  is  expressed 
by  formula  (11)  and  heat  balance  equation  (2)  acquires  the  following 
form: 


Designating 


Figure  3*  Graph  for  Determining  the  Thawing  Rate  of  the 
Frozen  Ground  Cylinder  When  the  Column's  Heat  Absorption 
is  Constant. 


In  order  to  make  it  easier  to  integrate  equation  (19)  we  replace 
the  function  /r  In 


with  the  simpler  dependency 


where  c and  f are  constant  coefficients 


By  substituting  (28)  into  (22),  we  find 


> L ;V  .J 

t = ^ ~ir'  ~'’«r'~rJ  + 


21  'c  f)  i7T7Arth  J/^TTry''  +c"; 


A.-** i. 


.V,  = ~ 


By  determining  the  integration  constant  from  the  initial  condition 
that  when  t = we  have  r = r^,  we  obtain  the  final  computation  equation 

for  calculating  the  growth  of  the  cylinder  when  the  brine  temperature 
is  constant: 

- - - r -I  „T  t\  2.V,  / T lA 

I (3D 

2 , - - - - ’ rv7±r(Ar  ">  1 ' x^t  Ar  'T  '<)  ] 

If  the  influence  of  the  changing  freezing  conditions,  for  instance 
an  increase  in  the  temperature  of  the  brine  or  the  filtration  flow, 
reduces  (melts)  the  frozen  ground  cylinder,  then  inequality 

k 

— < 0 occurs,  and  then  instead  of  (27)  we  obtain: 


V T"  ,-i»  O 

y ------- /•■  — i he  -A-i — -- 


■?[*•«»  i '1  - /! 


Therefore  when  the  cylinder  melts,  it  is  necessary  to  substitute 
expression  (32)  into  equation  (22)  instead  of  equation  (27);  then 
after  making  the  appropriate  transformations  and  determining  the 
integration  constant,  we  find  the  following  computation  equation  for 
the  speed  with  which  the  cylinder  melts  when  the  brine  temperature 
is  constant: 

r i ' 2 .1  h 

B,  -«(r;  -/•;  } - v (c  'r  ; J 
, • ' (3 

-^r,- Arctl.J  ^Sjr,) 


-10- 


n /■ 


From  formulas  (31)  and  (33)  it  follows  that  the  freezing  rate 
is  the  function  of  two  parameters:  B^  and  N^.  The  coefficients 

f and  c are  constant  and  are  determined  from  expression  (20);  thus, 
for  instance,  if  the  radius  of  the  column  is  r^  = 0.05  meters,  then 

we  can  assume  (up  to  the  radius  of  the  frozen  ground  cylinder  r < 0.5  m) 

-1/2 

that  c = 3.8  m and  the  coefficient  is  f = cr^  = -3. 8-0. 05  = 

1/2 

= 0.19  m . Dependencies  (31)  and  (33)  are  used  to  make  computations 
for  the  case  where  r^  = 0.05  m.  The  results  of  these  calculations 

are  shown  in  Figures  4 and  5,  and  they  make  it  possible  to  find  a 
solution  for  them  fairly  quickly  instead  of  finding  them  from  the 
formulas . 


Figure  4.  Graph  to  Determine  the  Speed  with  which  a Frozen  Ground 
Cylinder  is  Formed  at  Constant  Brine  Temperature. 


The  maximum  value  of  the  frozen  ground  cylinder  radius  can  be 
determined  if  the  left-hand  portion  of  equation  (17)  is  made  equal  to 
zero: 


r In 
max 


max 


= N. 


(34) 


-11- 


is 


Figure  5.  Graph  for  Determining  the  Speed  with  Which  the 
Frozen  Ground  Cylinder  Melts  When  the  Brine  Temperature  is 
Constant. 

The  results  of  calculating  r = f(N.J  for  r.  = 0.05-0.075-0.1  m 

& max  2 0 

are  shown  in  Figure  6. 


Figure  6.  Graph  for  Determining  the  Maximum  Radius  of  the 
Frozen  Ground  Cylinder. 

Thus  it  is  found  that  the  maximum  radius  of  the  frozen  ground 
cylinder  depends  only  on  parameter  N?  when  the  radius  of  the  column  is 

given.  In  order  to  analyze  the  influence  of  the  individual  factors 
on  the  maximum  cylinder  radius,  we  will  determine  the  value  of  N?. 


For  this  purpose  we  insert  expressions  (7)  and  (18)  into  (30), 
and  then  we  find 


N, 


1,39 


/ .i 


(:>  - iJ«i  ' 


(35) 


An  examination  of  this  equation  and  Figure  6 makes  it  possible  to 
draw  the  following  conclusions: 


1)  changes  in  the  brine  temperature  and  the  filtration  flow 
temperature  influence  the  course  of  freezing  more  significantly  than 
changes  in  the  flow  rate  of  the  filtration  flow; 


2)  when  water-permeable  ground  freezes,  the  maximum  cylinder  radius 
can  only  be  increased  by  reducing  the  brine  temperature  or  by  increasing 
the  column  radius.  The  other  values  are  usually  given  and  cannot  be 
changed; 

3)  the  maximum  radius  of  the  frozen  ground  cylinder  is  not  a 
function  of  the  amount  of  latent  ice  formation  heat. 


4.  Evaluation  of  Previous  Assumptions 


The  assumptions  adopted  in  this  article  have  been  previously 
introduced  by  many  authors  to  one  extent  or  another  [1,  2,  3,  4,  5], 
but  the  influence  of  these  assumptions,  as  far  as  we  know,  have  not 
been  evaluated  in  detail.  Therefore  the  evaluation  made  below  of 
these  assumptions  is  one  of  the  first  attempts  undertaken  in  this 
direction. 


First  assumption:  the  thermal  flux  which  changes  the  temperature 

of  the  frozen  ground  can  be  taken  into  account  by  increasing  the 
calculated  value  of  the  latent  ice  formation  heat. 


The  thermal  flux 


is  equal  to: 


(36) 


where  w is  the  heat  content  of  the  frozen  ground. 

The  heat  content  of  the  elementary  ring  of  the  frozen  ground 
cylinder  is  equal  to 


dm  — 2~c. -/, ) p>  df>, 

where 

c^  is  the  thermal  capacity  of  the  frozen  ground; 


(37) 


tp  is  the  temperature  of  the  elementary  ring  of  the  frozen  ground; 
is  the  percentage  of  frozen  ground; 


p is  the  radius  of  the  elementary  ring. 


As  is  known  from  heat  transfer  courses,  the  temperature  distribution 
established  in  the  cylinder  is  determined  by  the  following  equation  [6] : 


= S-lii-'-. 

f ' rn 

In  — 
ra 


The  temperature  of  the  brine  is  equal  to: 


0 — In  - - ■ ( 

2 -'.J.  r- 

By  substituting  (39)  into  (38)  we  find: 

f 0 2 t>.L  y 

then  for  the  heat  content  of  the  elementary  ring  of  the  cylinder  we 
find 

d~sS--  — !n  --  do. 

By  integrating  the  latter  expression  within  the  limits  of  o = r( 
to  p = r,  we  find  the  heat  content  of  the  frozen  ground  cylinder: 


W“  ' 4 " ' 

from  which  we  find 


V'f  - ;§-c Mi± 


If  r > 4r^,  then  it  is  obvious  that 


1""  100-i=G;.',: 


and  therefore  by  making  a minor  error  on  the  side  of  increasing  the  thermal 
flux  Q , we  assume  ['<  _ ^ and  therefore  we  will  have 

8 r ’ 


Q - - ir  . 
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This  expression  makes  it  possible  to  represent  the  sum  of  the  thermal 
fluxes  Qi  [expression  (3)]  and  Qg  [equation  (43)]  in  the  following 

form: 


tir 


QfH-Qg~(=-. 


where 


'8 


Azi.L 


(44) 


Designating 


we  finally  obtain: 


a'  = a + a , 
g 


Qf-Qs 


(45) 


(46) 


Expression  (46)  makes  it  possible  to  allow  for  the  thermal  flux 
which  reduces  the  temperature  of  the  frozen  ground  by  increasing 
the  calculated  value  of  the  latent  ice  formation  heat.  We  should 
emphasize  the  fact  that  the  freezing  time  of  the  ground  is  directly 
proportional  to  the  value  of  the  latent  ice  formation  heat,  but  the 
latter  exerts  no  influence  on  the  maximum  size  of  the  frozen  ground 
cylinder. 

An  examination  of  expressions  (44)  and  (45)  makes  it  clear  that 
for  practical  purposes  the  increase  in  the  calculated  value  of  the 
latent  ice  formation  heat  does  not  exceed  10-20%.  Consequently, 
the  value  in  (44)  can  be  determined  approximately,  and  this  has 

little  effect  on  the  final  results  of  calculating  freezing  time. 

Second  assumption:  at  every  moment  in  time  the  temperature 

distribution  in  the  frozen  ground  corresponds  to  the  established 
thermal  state. 

This  assumption  is  reflected  in  the  calculated  values  of  the 
column's  heat  absorption.  The  heat  absorption  of  the  column  is 

equal  to  the  product  of  three  values,  i.e.:  the  coefficient  of  the 
frozen  ground's  thermal  conductivity,  the  area  of  the  column's  lateral 
surface  and  the  temperature  gradient  in  the  frozen  ground  near  the 
column  surface. 

1 

The  temperature  gradient  in  the  ground  near  the  column  wall  is 
always  less  than  when  the  thermal  state  is  established  than  when  it  is 
not  (Figure  7).  Therefore  the  above-mentioned  assumption  about  the 
established  temperature  distribution  in  the  frozen  ground  should 
mean  that  the  calculated  heat  absorption  values  of  the  column  will  be 
below  the  actual  values.  Consequently,  the  heat  balance  equation 

(1)  Qf  = Q - Qf  * Q provides  understated  values  of  the  thermal  flux 

8 
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Figure  7-  Temperature  Distribution  in  the  Cylinder. 

1,  Established  temperature  distribution;  2,  unestablished 
temperature  distribution. 


However,  it  is  obvious  that  the  difference  between  the  actual  value 
of  the  column's  heat  absorption  and  the  calculated  value  is  less  than  the 
value  of  the  thermal  flux  which  changes  the  temperature  of  the  frozen 

ground.  This  allows  us  to  state  that  if  we  use  heat  balance  equation 
(2)  Qf  = Qe  - Qj-.  then,  quite  the  contrary,  understated  calculated  values 

will  be  obtained  for  the  thermal  flux  from  the  latent  ice  formation  heat, 
and  consequently  understated  calculated  values  for  the  freezing  rate. 

The  true  value  lies  between  the  two  calculated  values. 

The  relative  difference  between  the  two  solutions  is  equal  to: 


j -l  r 


g 


(47) 


and  because  of  what  is  stated  above  concerning  the  first  assumption, 

this  difference  is  found  to  be  within  20%  for  practical 

purposes. 

Third  assumption:  frozen  ground  has  the  shape  of  a cylinder. 

The  reason  for  the  circular  shape  to  be  violated  may  be  non-uniform 
distribution  of  the  thermal  influx  from  the  unfrozen  ground  over  the 
cylinder  surface. 

The  non-uniform  distribution  of  the  thermal  influx  is  characterized 
by  the  relationship: 


= ^fi 


^fi  av 


(48 ) 
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) 


I 


I 


i .. 


where 


is  the  intensity  of  the  thermal  influx  from  the  unfrozen  ground 
at  a given  point  on  the  cylinder  surface; 

q_.  is  the  average  intensity  of  the  thermal  influx  from  the  unfrozen 

I X XV 

ground . 

We  will  determine  the  intensity  of  the  thermal  influx  (the  position 

of  the  coordinate  axes  relative  to  the  direction  of  the  filtration  flow 
is  shown  in  Figure  8). 


Figure  8.  Frozen  Ground  Cylinder  in  the  Filtration 
Flow. 

By  using  the  solution  to  the  task  of  the  cooling  of  the  laminar 
flow  when  the  flow  around  the  cylinder  is  transverse  [7],  with  regard 
to  our  case,  we  can  write 


‘/fi 


' ifW  w I-  W — !>o) 
r ifr,  ’ 


where  is  the  current  function  of  the  filtration  flow. 
The  current  function  n is  expressed  by  the  equation 

r,*=vl2r—x  -'-—V 

X1  ■■■  V-  ! 


(49) 


On  the  cylinder  surface  where  ]/  xl  -j-y-  = r,  we  have: 


or,  noting  that  x/r  = cos  a where  a is  the  radius-vector  argument  of  the 
point  on  the  cylinder  surface  (Figure  8),  we  find 
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r,  = 'lvr  (l  --cos  a). 


(50) 


By  substituting  expression  (50)  into  equation  (49)  and  integrating 
in  terms  of  n»  we  obtain: 

Q = j ?firfr(  = 2|/  (ft— »„)  /-  \/T~cZT,  -- 

where  is  the  integration  constant. 

By  differentiating  the  latter  expression  in  terms  of  a and  by 
substituting  the  result  into  equation 


0 _ j-Q_ 
'“ft  Lrch’ 

we  will  have:  , — 


or,  after  certain  trigonometric  transformations,  we  will  finally  obtain 


?fi  = 2] 


(51) 


The  average  intensity  of  the  thermal  influx  from  the  unfrozen  ground 
is  equal  to: 

flU  — ?-L. 

ft  av  U-rL 

Substituting  expression  (4)  into  this  expression,  we  obtain: 


%av~  — }/  ~Tr~  (•> — ”o>- 
By  substituting  (51)  and  (52)  into  (48),  we  find: 


(52) 


2l  = — cos  — 
2 2 


(53) 


The  dependency  = f(a)  is  plotted  in  Figure  9.  An  examination  of 

it  shows  that  the  thermal  influx  from  the  unfrozen  ground  is  distributed 
along  the  surface  of  the  cylinder  in  an  extremely  uneven  fashion.  This 
fact,  however,  does  not  violate  the  circular  form  of  the  cylinder  since 
the  half-sum  of  the  thermal  influx  intensities  at  diametrically  opposed 
points  on  the  cylinder  surface  is  equal  to  the  average  intensity  of 
the  thermal  influx. 
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This  condition  is  written  in  the  following  fashion: 


The  dependency  = f(a)  is  shown  in  Figure  9.  An  examination  of 

it  shows  that  the  values  of  deviate  from  unity  by  no  more  than  22%. 

From  this  it  follows  that,  despite  the  extremely  uneven  distribution 
of  the  thermal  influx  from  the  water-permeable  ground,  the  frozen 
ground  has  a shape  which  is  similar  to  that  of  a cylinder,  but  its 
axis  is  shifted  in  the  direction  of  the  flow:  an  eccentricity  arises 

between  the  cylinder  axis  and  the  axis  of  the  freezing  column. 


Figure  9-  Distribution  of  Heat  Influx  from  Unfrozen  Ground 

Over  C i rcumference  of  Cylinder. 

The  correctness  of  this  conclusion  has  also  been  confirmed  by  the 
carried  out  by  the  Ice  Heat  Laboratory  of  the  B.  Ye.  Vedeneyev  All-Union 
Scientific  Research  Institute  of  Hydraulic  Engineering. 

Figure  10  shows  a photograph  of  an  frozen  ground  cylinder  formed 
in  sandy  ground  at  a flow  speed  of  v = 0.06  m/hr  and  a water  temperature 
of  d = 10. 6°C. 


Figure  10.  Frozen  Ground  Cylinder  (the  Filtration  Flow 
was  Directed  From  Right  to  Left). 

The  fourth  assumption:  the  axis  of  the  frozen  ground  cylinder 

coincides  with  the  axis  of  the  freezing  column. 

The  essence  of  this  assumption  consists  of  the  fact  that  in 
calculating  the  heat  absorption  of  the  column  it  is  possible  to 
ignore  the  eccentricity  between  the  axes  of  the  cylinder  and  the 
freezing  column.  In  order  to  determine  the  error  which  is  committed 
in  this  process,  let  us  investigate  the  relationship 


(gc),-„ 


(56) 


where  e is  the  eccentricity  between  the  cylinder  and  the  freezing 
column. 


As  is  well  known,  when  there  is  an  eccentricity  [8] , 


-from  this  we  find 


(57) 


(58) 


r7 


i 


From  Figure  11,  which  is  plotted  with  the  aid  of  dependency  (58) , 
it  follows  that  since  e < 0.5  r,  then  the  error  which  is  committed 
in  calculating  the  thermal  absorption  of  the  column  without  allowance 
for  the  eccentricity  is  small  (less  than  10%). 

Thus,  we  may  consider  that  it  has  been  proven  that  the  introduction 
of  these  four  assumptions  introduces  no  large  errors  into  the 
calculations  of  the  freezing  of  water-permeable  ground. 


Figure  1 1 . 


Function 


5.  Sample  Calculation 

Let  us  show  by  way  of  an  example  how  we  use  the  calculated  dependencies 
obtained  here. 

The  task  may  be  formulated  as  the  following:  it  is  necessary  to  form 
a frozen  ground  cylinder  of  radius  r = 0.7  meters  in  a filtering  layer 
of  ground  10  m thick.  The  cylinder  is  set  perpendicular  to  the  direction 
of  the  flow. 

The  following  data  are  known: 

speed  of  filtration  flow 

temperature  of  filtration  flow., 
freezing  temperature  of  ground., 
thermal  conductivity  factor  of 
unfrozen  ground 

thermal  conductivity  factor  of 
frozen  ground 

latent  heat  of  ice  formation  per 

unit  of  ground  volume 

thermal  capacity  of  water 

specific  weight  of  water 

radius  of  freezing  column 

length  of  freezing  column 


v = 0. 1 m/hr 

0 = 3°C 

d = - 1°C 
0 

= 1 kcal/m-deg'hr 

= 2 kcal/m-deg-hr 

o = 24,000  kcal/m-5 
c^  = 1 kcal/kg-deg 

Yw  = 1,000  kg/mJ 

Tq  = 0.10  m 

L = 10  m. 
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The  cooling  device  and  its  operating  mode  are  selected  in  such 
a way  that  the  heat  absorption  of  the  column  changes  as  the  size  of  the 
frozen  ground  cylinder  increases,  as  shown  in  Figure  12;  in  this  process 
the  temperature  of  the  brine  does  not  reach  a constant  maximum  low 
temperature. 

The  following  have  to  be  determined: 

1)  the  time  required  for  the  formation  of  a frozen  ground  cylinder 
of  the  given  dimensions; 

2)  the  column's  heat  absorption  value  required  to  maintain  the 
frozen  ground  cylinder  in  the  frozen  condition  (passive  freezing); 

3)  the  brine  temperature  in  the  column  during  passive  freezing. 


Ilil'M 


5 S’  C~  C-3  03  1.0  «- 

Figure  12.  Change  in  Heat  Absorption  of  Freezing  Column. 


Solution.  We  will  determine  the  freezing  time.  First  according  to 
expression  (7)  we  find  the  value  of  parameter  A: 


8|/L 


1-1000.0,1 


■WP+  1)-= 


= 1810  kcal/hr-m 

By  substituting  the  initial  data  into  formula  (8),  we  find  that  the 
freezing  time  is  determined  by  the  following  expression: 


~ It.Lz  10-24  000  l‘ - = dr  - 

.)  i(r)-AVr  J r 


1.51  •10s  1 ' — dr. 

J / (r)  - 1810  J r 


We  will  carry  out  the  calculation  by  the  tabular  integration  method, 
for  which  purpose  we  divide  the  integration  limits  into  n sections:  each 
section  will  be  designated  by  a serial  number  i.  Then  we  will  have 
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r,  = 1,51-100  V =-  Ui-r,,). 

i 1 

The  value  of  the  column's  heat  absorption  will  be  selected  according 
to  Figure  12,  allowing  for  the  fact  that 

The  rest  of  the  calculations  require  no  explanation  and  are  summarized 
in  the  following  table: 


/ 

JZ 

N 

3 " 

V. 

0 '■N 

M 

- JZ 

o ; 

E j* 
— 

j£  | 

© ^ 

■ ; 

o 1 

* -X  \ 
< 1 

73 

2 2 
^JZ 

i * 

*«S  3 

; 

at 

_ * _ _ 

l<>73 

Jr- 

JO 

r 

i 

1 •!-.  hr 

>— 

I 

0,15 

720J 

700 

0500 

23,08-10-'* 

1 

i 0,05 

1,15-10-® 

2 

2 

2 

0,25 

5900  1 

905 

4995 

50-10- u 

1 0,10 

3-10-“ 

8 

10 

3 

0,35 

4930 

1070 

! 3880 

90,3-10-" 

i 0,10 

9,03-10-® 

14 

24 

1 

0,45 

1400 

1210 

1 3190 

141  • 10~® 

0,10 

14,1-10® 

21 

45 

f, 

0,55 

1100 

: 1.340 

2700 

19S-10-® 

0,10 

j 19,8-10  ® 

30 

75 

6 

0,70 

37.70 

, 

j 1510 

2210  ; 

1 

.312-10-® 

0,15 

1 46,8-10-® 

71 

110 

Commas  indicate  decimal  points. 

The  calculation  results  show  that  the  formation  of  an  frozen  ground 
cylinder  wiil  take  146  hours,  i.e.,  approximately  six  days  (Figure  13). 

Let  us  compare  this  result  with  that  which  is  obtained  if  we  carry 
out  the  calculation  according  to  formula  (15). 

The  column's  heat  absorption  drops  from  an  additional  value  of 
8,000  kcal/hr  to  3,600  kcal/hr  at  the  beginning  of  the  freezing  process 
(see  Figure  12).  We  will  assume  that  the  column's  heat  absorption  is 
a constant  value,  equal  to 

SOOO :J6oo 

Qc = 5800  kcal/hr. 

From  expressions  (13)  and  (14)  we  find: 


„ 2-t.i  2-- 10-21  non 

B,  — ( nr  ...  . 

A ISIO 

v _ !Ls . i8_°L>  „ n o 

‘ ‘ A 1810  "'■*  ■“ 

Then,  according  to  formula  (15),  we  obtain 


- 0,1  -)  4-  3,2  (r,  - 0,1)  -f-  2-3, 2-{r.r 
± 1 

3,2  — r,~ 


2- 3,23  In- 


3,2  - O.t  • 


0,1  2/-f 


for  ^ = 0.7  m we  find  t = 84  hours. 

The  results  obtained  previously  (by  the  tabular  integration  method) 
t = 146  hours  is  more  precise.  Thus,  we  see  that  averaging  the  column's 
heat  flow  within  such  large  limits  and  for  the  entire  freezing  period 
provides  only  a very  approximate  result  (r  = 84  hours).  It  is  much 
more  precise  to  determine  the  freezing  time  in  sections,  assuming  that 
each  section  has  its  own  averaged,  constant  value  for  the  column's 
heat  absorption.  We  will  demonstrate  this  by  means  of  the  same  example. 


In  carrying  out  the  calculations 
we  will  utilize  the  graph  shown  in  Figure 
2. 


We  will  divide  the  desired  freezing 
time  into  four  parts: 


• — ri  ~ "»  ' "a  t.,. 

1/2 

As  before,  A = 1810  kcal/hr-m  and 
Bx  = 834  hr/m3/2. 

We  will  first  determine  the  formation 
time  (Tj)  for  a frozen  ground  cylinder 

0.3  m in  radius. 

The  average  heat  absorption  value 
is  equal  to  (Figure  12) : 

_ M00  + 4WO  kcal/hr. 

c 2 


Figure  13-  Graph  of  Freezing 
Speed . 


and  consequently 


N«  = 


£c 

A 


1,100 

1810 


- 0,0  .x 


1/2 

According  to  Figure  2 we  find  that  when  r = 0.3  m and  = 3.5  m , we 
have 


L « 0,01  .11  - . 

B. 


from  which  we  obtain 


*i  — 0,'MC,  «■  t.o,. il  % . hrs. 


Now  we  will  determine  the  freezing  time  from  r^  = 0.3  m to 


r2  = 0.5  m. 


The  average  heat  absorption  of  the  column  is  equal  to: 


4 SCO  4-  4W) 


4ii'0  kcal/hr. 


and  then 


1S10 


From  Figure  2 we  find  that  when  r = 0.5  m and  = 2.43  m , we  have 


— - 0,07  .«  ' . 


From  the  same  Figure  2,  we  find  that  when  r = 0.3  m and  - 2.43  m ^ we  have 


Then  we  find 


0.02  .11 * . 


M4  (0  0"  — 0,02)  42  his. 


For  the  third  freezing  section  from  r^  = 0.5  m to  r^  = 0.6  m. 


we  find: 


0 w • 5000  kcal/hr. 

• c 1 


.V,  ~ — - — - 2.1.)  .«  - . 
lbl'> 


From  Figure  2 we  determine: 


for  r. 


• 0.12  If  ' . 


and  for  r. 


--  - 0.0s  .n  - 
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I 

* 


Consequently, 


*3  - 334  (0,12  — 0,03)  - 33  hrs. 


Let  us  determine  the  freezing  time  from  r^  = 0.6  m to  r^  = 0.7 
We  find: 


_ 3300  - 3000.  ,j37OT  kcal/hr 
9 


,,  3roo  - 

A , — =■  2,0o  m 

1310 


According  to  Figure  2 we  have: 


for  r„ 


— = 0,19  .»  ' 


for  r, 


Consequently 


--  = 0,13  ,ii  ' . 
B i 


_ s-u  (0,19  - 0,13)  - SO  hrs. 


Finally  we  find  that  the  freezing  time  is  equal  to 

- - 3 - 42  - 33  - 30  = 133  hrs. 

This  result  is  considerably  more  precise  than  the  previous  one 
(t  = 84  hours)  since  it  makes  sufficient  allowance  for  the  change  in  the 
column's  thermal  absorption  during  freezing. 

Let  us  solve  the  second  portion  of  the  task:  determining  the  value 
of  the  column's  heat  absorption  during  passive  freezing. 

The  following  equality  should  be  valid 

Qc  = V 

According  to  expressions  (4)  and  (7), 

V’-1' 


In  our  case  A = 1810  kcal/hr-m  and  r = 0.7  m.  We  find  that  in  order 
to  keep  the  frozen  ground  cylinder  in  a frozen  state,  the  heat  absorption  of 
the  column  must  be: 

-26- 


Qc„  .4  I ~r  -=  16*0  1 0.7  =»  1510  kcal/hr. 


In  order  to  answer  the  third  question  of  what  the  temperature  of 
the  brine  during  passive  freezing  should  be  equal  to,  we  will  utilize 
formula  (11),  according  to  which 


m 


By  substituting  our  known  data  into  this  expression,  we  find 


'nr 


1510  , 0.7 

— In  -- 

2--M0  0,1 


— 2t.rc. 
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